Phase-shift inversion in oscillator systems with periodically switching couplings 
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A system's response to external periodic changes can provide crucial information about its dynam- 
ical properties. We investigate the synchronization transition, an archetypical example of a dynamic 
phase transition, in the framework of such a temporal response. The Kuramoto model under pe- 
riodically switching interactions has the same type of phase transition as the original mean-field 
model. Furthermore, we see that the signature of the synchronization transition appears in the 
relative delay of the order parameter with respect to the phase of oscillating interactions as well. 
Specifically, the phase shift becomes significantly larger as the system gets closer to the phase tran- 
sition so that the order parameter at the minimum interaction density can even be larger than that 
at the maximum interaction density, counterintuitively. We argue that this phase-shift inversion is 
caused by the diverging relaxation time, in a similar way to the resonance near the critical point in 
the kinetic Ising model. Our result, based on exhaustive simulations on globally coupled systems as 
well as scale-free networks, shows that an oscillator system's phase transition can be manifested in 
the temporal response to the topological dynamics of the underlying connection structure. 

PACS numbers: 05.45.Tp, 05.45.Xt, 64.60.an 



I. INTRODUCTION 

A collection of interacting coupled oscillators is one of 
the most intensively studied systems showing a dynamic 
phase transition called a synchronization transition [H-Q . 
Such systems have been used as models to describe vari- 
ous phenomena such as the pacemakers of the heart, the 
collection of amoeba, large-scale ecosystems and 
fiber-optic networks of optoelectronic oscillators l3| • As 
the coupling strength increases, such systems undergo 
a phase transition from the desynchronizcd (disordered) 
state to the synchronized (ordered) state characterized 
by critical phenomena. The most well-known theoretical 
formalism is the Kuramoto model [l[ where the oscil- 
lators interact with each other by a sinusoidal coupling 
strength with respect to the phase differences of the two 
oscillators. 

Aside from the steady-state behavior after the initial 
transient behavior, the dynamical aspect of the oscil- 
lators has also been studied recently [ll], [T^]. Espe- 
cially, the effects of temporally varying interaction struc- 
tures themselves are worthwhile to investigate since 
we can systematically analyze the response of sy stems to 
such structural changes of interactions [3, [isf . For in- 
stance, the mobile oscillator is considered in Refs. [T6l - [l8j 
as an example of temporally switching interactions. In 
this Brief Report, we take the Kuramoto model with pe- 
riodically switching interactions and analyze the tempo- 
ral response in the collective phase of oscillators to such 
periodicity in interactions. Such periodically switching 



interactions may have implications to biological or social 
communications among compartments under natural cir- 
cadian rhythms, for instance [l9| . 

In terms of the temporally averaged order parame- 
ter and its amplitude of oscillation, we veriiy that the 
model undergoes the same synchronization transition as 
the mean-field (MF) Kuramoto model [3[ with the finite- 
size scaling (FSS) analysis. Furthermore, with numerical 
simulations and analogy to the kinetic Ising model j20| . 
we show that the signature of the synchronization tran- 
sition is revealed in the relative phase shift of the order 
parameter with respect to the density of interactions, as 
a resonance-like phase-shift inversion. The result is also 
consistent with the diverging relaxation time at the crit- 
ical coupling strength [TllllJl- Therefore, we claim that 
the system's temporal response to the temporally chang- 
ing interactions can be used as an indicator of a phase 
transition. 



II. MODEL 

We consider the Kuramoto-type oscillator dynamics Q 
composed of N oscillators as 



(1) 
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where is the set of oscillators' phases, {u!i} is the set 
of natural frequencies of oscillators given by the Gaus- 
sian distribution with average and unit variance, and 
the symmetric adjacency matrix Oy- = 1 if z and j are 
connected and otherwise. The interaction a^Ty (t) me- 
diated by the edge between the oscillators i and j is sub- 
ject to the periodic function as (t) = H [sm{nt + 9ij )] , 



where H is the Heaviside step function, n is the frequency 
of switching common to all the edges, and {0ij} for i~j 
pairs with = 1 is the set of edges' intrinsic phases, 
which is randomly assigned from the uniform distribu- 
tion of the interval [0,7r). By choosing such an interval, 
the edge densities change periodically with the exact tri- 
angle wave subject to the frequency f2, from fully con- 
nected edges to isolated oscillators without interactions. 
Note that the factor 2 in Eq. ^ is to directly compare 
the coupling strength K of our model with T^j = 1/2 to 
the K value in the original Kuramoto model, where x 
refers to the temporal average of x(t) over a period. 

In addition to the original control parameter of cou- 
pling strength K in the Kuramoto model, the edge fre- 
quency ri driving the system is a crucial parameter as 
well. The response of oscillators to the toggling interac- 
tions is measured by the temporal behavior of conven- 
tional phase order parameter A(t) where 



1 ^ 

A(t) exp[i^(<)] - - ^cxp[i0,(t)] 



(2) 



i=i 



The order parameter's dynamics depends on the pe- 
riodic change of interactions with the density p(t) = 
J2t<j'^tjTij{t)/Yl,i<jaij. We provide the Java applet 
simulating this model system in case of globally con- 
nected oscillators, for readers to observe the behavior of 
a system of relatively small sizes |2l| . From now on, we 
present our numerical simulation results and their impli- 
cations. 
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p(t > 20) 

FIG. 1: A typical belmvior of the order parameter A(t) and 
the interaction density p{t). The system size A'' = 800, the 
edge frequency Q = 1, and the time series of A{t) depending 
on K are shown as smooth curves in (a) and p{t) (independent 
of K) as black triangle waves in (a). The A-p diagram for 
t > 20 after the transient behavior is shown in (b). All the 
curves are results averaged over 50 samples. 



III. RESULTS 
A. Globally coupled oscillators 

First, we consider the globally coupled case, i.e., Uij = 
1 — Sij . A typical temporal behavior of oscillators is il- 
lustrated in Fig. [TJ in case oi ft = 1. Note that the 
oscillating behavior of the order parameter A(t) driven 
by switching edges (with the same frequency) is shown 
even in case of quite low values of K , for which the orig- 
inal MF Kuramoto model corresponds to the disordered 
phase. Another interesting aspect is the relative phase 
shift of A{t) with respect to p{t). In other words, the 
maximum or minimum of A{t) occurs after some time 
when the maximum of minimum of p{t) is reached, as 
shown in Fig.[TJa). Moreover, as K is increased, the delay 
for maximum values and the delay for minimum values 
become more asymmetric, as clearly shown in Fig.[TJb). 
In spite of such temporal oscillations of A(t). the tempo- 
rally averaged value A (averaged after the transient time 
~ 20) shows exactly the same universality class of syn- 
chronization transition as the MF Kuramoto model with 
13 = 1/2, and v = 5/_2 as described in the FSS scal- 

ing collapse with A = N-^^" f {{K ~ K^in)) N^^^;^) 
(the scaling function / itself, as well as the critical cou- 
pling strength Kc, depends on 17) shown in Fig. We 
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FIG. 2: FSS scaling collapse of the temporally averaged order 
parameter A for various system sizes with the critical expo- 
nents of the original MF Kuramoto model [Kc = 1.62(1), 
/? = 1/2, and i/ = 5/2]. The edge frequency is = 1. 

suggest the slightly larger value of ifc(f^) than that of 
the MF model K^^ = ~ 1.596 stem from ad- 

ditional quenched randomness caused by {Oij} for edges, 
similarly to the lower Tc (the disordered phase extended) 
for the kinetic Ising model under the external oscillating 
field [lOl. Since we have confirmed the MF transition 
for A, from now on we focus the temporal behavior of 
A(t), especially with respect to the oscillating interac- 
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FIG. 3; The difference A between maximum and minimum 
values of A{t) oscillation and the area S enclosed by the loop 
in A-p the diagram are shown in the K-Q plane. The system 
size A*' = 800, and the color-coded values are A (a), A^l (b), 
S (c), and SO. (d). The time series A(t) for t > 20 after the 
transient behavior is used, and all the results are averaged 
over 50 samples. 
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FIG. 4: AQ (a) and SQ (b) as functions of K, for various 
values. The black vertical lines correspond to the MF critical 



coupling strength 



= 2, 



tion strength p{t). 

To characterize the oscillating behavior more quantita- 
tively, we systematically measure the difference between 
maximum and minimum values of A(t) oscillation de- 
noted as A (roughly twice the amplitude of A(i)) and 
the enclosed area in the A-p diagram denoted as S, for 
a wide range of $7 and K values. Note that the data 
shown here are from the system size = 800, but we 
have cheeked that there exists no significant finite-size 
effects based on the simulations results with smaller and 
larger systems. The time series A{t) for t > 20 after the 
transient behavior is used, and all the results are from av- 
eraged over 50 samples for all the simulations from now 
on. As shown in Fig. [3l the scaling behavior of A and 
S is quite similar, and notably both A and S seem to 
be inversely proportional to ft for O > 1.2. This scaling 
behavior makes the and Sfl in Figs.[31[b) and (d) de- 
pend only on K ioi ^ > 1.2, and is also clearly observable 
from the fact that the curves for > 1.2 are collapsed in 
Fig.m Since f2 represents the angular velocity of the fiuc- 
tuation of edges (hence that of A(t)), Aft corresponds to 



the "linear velocity" of A(t), which is shown to be con- 
served for a given value of K. More importantly, the 
synchronization transition is clearly shown in the steep 
change of A and S near the critical coupling strength 
for the original MF Kuramoto model, as shown in 
Figs. [3] and 131 Therefore, we conclude that the signature 
of the MF synchronization transition of Kuramoto model 
is resurfaced in terms of A and 5*. 

Aside from the amplitude of A and enclosed area in the 
A-p diagram, a notable property of the A-p diagram in 
Fig. [TJb) is successive change in the shapes of enclosed 
regions, as K is increased. First, the desynchronization 
process of A{t) after the onset of decreasing phase of p{t) 
is slower than the synchronization process of A{t) after 
the onset of increasing phase of p{t) in general. This 
asymmetry is natural, since the response of the order 
parameter to such linear functions is integrated form, i.e., 
quadratic functions depending on the initial values. 

We show that the synchronization transition is also 
observed in the phase shift of A{t) with respect to p{t) 
as well. To systematically analyze such a phase shift 
caused by the temporal delay of A(i), we check the rel- 
ative height difference of A at the right and left ends 
in Fig. [Hb), corresponding to A(/9max) and A{pmin) re- 
spectively, also affects the shape of enclosed regions. In 
Fig. [SJa), we plot such height differences in the K-ft 
space, and the height difference occurs mainly in small 
57 and large K regimes. Interestingly, near K ~ 1.6 in 
relatively large f2 regimes, the difference becomes nega- 
tive, which corresponds to A(pmax) < A(pniin)- We call 
this phenomenon phase-shift inversion. .ftT ~ 1.6 is close 
to the MF critical coupling constant again, and we 

suggest that it is related to critical slowing down near 
the critical point in a phase transitiorr, as shown in the 
diverging relaxation time r at [Tl|, [T^], given by 

1 = y/^Kexp[l/{2T'^)]eik[l/{V2T)]. This resonance- 
like phase-shift inversion phenomenon occurs similarly 
to the frequency matching condition around the phase 
transition in the MF kinetic Ising model [l^]. In other 
words, the MF synchronization transition in the Ku- 
ramoto model is also reflected in the temporal response 
to the periodically switching interactions. 



B. Scale-Free Networks 

To validate our conclusion on the phase transition in 
terms of phase-shift inversion, we use substrate network 
topologies other than the globally coupled one. The 
static scale-free network (SFN) model [23] with the tun- 
able degree distribution p{k) ^ k~'^ provides various un- 
derlying interaction structures such as the pure MF case 
(7 > 5), the MF with critical exponents depend on 7 
(3 < 7 < 5), and eventually the one where critical cou- 
pling strength vanishes as i^c — J' (7 < 3) [2^. We 
use the SFN model for {uij} in Eq. ([Ij and observe the 
phase-shift inversion diagram for the three representative 
cases: 7 = 6.5, 3.75, and 2.2, as shown in Figs. [5Jb)-(d). 
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FIG. 5: Height difFerence between A(pmax) and A(/9inin), for 
the globally coupled case (a) and SFNs with 7 = 6.5 (b), 3.75 
(c), and 2.2 (d). Note that the rescaled coupling strength 
KN/{k), where (k) is set as 20, is used for (b)-(d), to compare 
with the globally coupled case (a) with the same Eq. The 
same system size A'^ = 800 is used for all the cases. 



same mean-field synchronization transition as the origi- 
nal Kuramoto model with the finite-size scaling analysis, 
and found that the phase transition is observed in terms 
of dynamical properties such as the amplitude of oscilla- 
tion of the order parameter and its relative phase shift 
with respect to the overall strength of interactions. In 
particular, the latter causes the phase-shift inversion phe- 
nomenon that the significantly large phase shift near the 
phase transition let the order parameter at the minimum 
interaction density to be larger than that at the maxi- 
mum interaction density. The conclusion does not only 
hold for the globally connected substrate structure, but 
also holds for other types of substrate structures such as 
scale-free networks. Such observations strongly suggest 
that our model is related to the mean-field kinetic Ising 
model with the frequency matching condition around the 
phase transition [20| and would be worth investigating 
further. Since periodicity plays important roles in many 
parts of the nature, especially for systems under natural 
circadian rhythms or external stimuli, our results suggest 
that the temporal response to such periodicity can be an 
important cue to characterize the system. 



One can observe the phase-shift inversion in SFN cases 
where the synchronization transition exists (7 > 3) as 
well, at lower coupling strength explained by Kc oc 1/ (fc^) 
when (k) is fixed [l^, where (x) represents the averaged 
X over the entire nodes in a network. The 7 < 3 case 
reveals its vanishing critical point in the fact that the 
seeming phase-shift inversion occurs in much lower values 
(and over the very narrow region) of coupling strength 
[Fig.EJd)], due to the finite-size effect. Therefore, the re- 
sults from SFNs also clearly support the fact that phase- 
shift inversion reflects the phase transition. 

IV. SUMMARY AND DISCUSSIONS 

We have studied the Kuramoto model under period- 
ically switching interactions, checked that it shows the 
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